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1. Introduction 



The role of the effective potential V(<p) in determining the nature of the vacuum in 
renormalisable field theories was emphasized in the classic paper of Coleman and Wein- 
berg[l] (CW). Their particular interest was in the special case when the renormalised value 
of the (mass) 2 parameter, m 2 , of the scalar fields was zero. They were able to demonstrate 
the occurrence of spontaneous symmetry breaking through radiative corrections for this 
case, as long as gauge fields are present. This scenario however, is excluded by current 
experimental limits on mu (the Higgs mass) and rrit (the top mass). 

In the m 2 < case, radiative corrections are still important in determining whether 
the tree minimum of V corresponds to the true ground state of the theory. ( For a review 
and references, see ref. [2] ). The one- loop Yukawa coupling contribution to V tends to 
destabilise the vacuum, and consequently leads to an upper bound for rat as a function of 
ra H [3]. 

The analysis described in ref. [3] was based on a renormalisation-group (RG) "im- 
proved" form [1,4] of V including one- loop corrections. Here we present the results of the 
two-loop radiative corrections to V, V^ 2 \ with a view to possible refinement of the bounds 
described in ref. [3] . 

Another motivation for our calculation is as follows: Given the standard model 2-loop 
/^-function, the dependence of V ^ on the renormalisation scale /i may be readily inferred 
from the fact that V satisfies a RG equation. (This calculation was in fact attempted in 
ref. [5]; we will comment later on this paper). We shall instead calculate the full V^ 2 \ using 
minimal subtraction (MS) throughout, and use the results to infer the 2-loop /3-functions 
for m 2 and the quartic Higgs coupling A ((3^). As we shall see, this procedure will 

(2) (2) 

expose some minor errors in the expression for the SM (3 X given in ref. [6]. As for (3 2 , as 
far as we know it has not previously appeared in the literature. As already mentioned in 
ref. [7], it is essential that be calculated using MS in order that be consistent with 
the RG equation with the MS /3-functions; the /3-functions in a multi-coupling constant 
theory are scheme dependent at the two-loop level [8]. Thus although there exists [9] (in 
the special case (h = m 2 = 0)) a calculation of V^ 2 \ this cannot be used for comparison 

(2) 

with the standard model p\ , since MS was not employed. 

The plan of this paper is as follows. In section 2 we describe briefly the effective 
potential formalism, give the result for the one-loop correction, V^ l \ and introduce some 
notation. In section 3 we explain our procedure for finding within the MS scheme, and 
in section 4 we apply a novel method to the evaluation of the pertinent Feynman integral. 
The method involves the use of differential equations [10] with the added refinement of the 
method of characteristics. Section 5 consists of a summary of the results, while in section 
6 we show how the fact that V(4>) satisfies a RG equation leads to a determination of j3^ 2 \ 

(2) 

and j3\ . Finally in section 7 we discuss applications. 
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2. The one-loop effective potential. 



The effective potential formalism of CW and the functional refinements introduced by 
Jackiwfll] are too well known to require any but the briefest review here. In general one 
shifts scalar fields as follows: 

0(a;)-0 + ^(x) (2.1) 

where is x-independent.Then the effective potential V(4>) is given by the sum of vacuum 
graphs, with ^-dependent propagators. Equivalently, one can calculate graphs with a single 
4> q external field, which, it is easy to show, leads to a determination of dV/d<p. The latter 
method is perhaps simpler at one loop, but not beyond since it leads to more Feynman 
diagrams. (Here we differ as a matter of opinion from ref. [9]). 

In the SM case, one can exploit gauge invariance to perform the shift of Eq. (2.1) on 
one only of the four scalar fields: 




G±{X) \ (2 2) 



We must also choose a gauge; the 'tHooft-Landau gauge is the most convenient one. (In 
fact V(4>) is gauge invariant only at its extrema; this gauge is a good one in the sense of 
ref. [12]). In this gauge the W, Z and 7 propagators are transverse, and the associated 
ghosts are massless and couple only to the gauge fields; the "would be Goldstone" bosons 
G^,G have a common mass deriving from the scalar potential only. 
At the tree level the effective potential is V^ ' (0), given by 

^ (O) (0) = ^0 2 + ^ 4 (2-3) 

(Note that our definition of A differs by a factor of 3 from that of ref. [6].) 
The result of the one loop calculation is: 

KV m = ^(Eff - §) + ?f!(EG - §) - 3T*(Er - §) 

+ 2^!(taw-5 ) + 5|!(Ez-5) (2.4) 

where 

k = 16tt 2 , H = m 2 + -c/ ) 2 ,T=— </> 2 

2 2 

G = m 2 + ^\W = ^ 2 , Z =^±A> 2 

— X 
and \nX = In — ^ + 7 — ln47r, 

7 being Euler's constant. (Note that the sign of the ln47r term is rendered incorrectly 
in ref. [7], Eq. (4)). Here h is the top quark Yukawa coupling (we neglect other Yukawa 
couplings throughout). 
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If m 2 < 0, then at the minima of V(°\<j)) we have G = and H, T, W, Z become 
the tree level (masses) 2 of the Higgs, top quark, W and Z bosons respectively. The 
non logarithmic terms in Eq. (2.4) may be altered (or indeed removed) by a change in 
subtraction scheme. Since, however, we wish to consider the RG equation with the MS 
scheme, it is essential that we retain them. ( The MS scheme corresponds, of course, to 
simply replacing \nX by ln(X/fx 2 ) throughout; this would not affect the RG analysis since 
all RG functions are identical in MS and M,S[13]). 



3. The two-loop calculation: preliminaries. 

In this section we outline our basic strategy for the calculation. By elementary ma- 
nipulation we can reduce each individual Feynman diagram to a sum of integrals of the 
form of either: 

T( \-^L[ ddkddq /o n 

{27v)d J {k 2 + x){q 2 + y){{k + q) 2 + z) ^.±j 

J(x,y) = J(x)J(y) (3.2) 



or 



J{x) = ^[J^ = (Sr T(1 * )x ^ (3 .3) 



where 

(2n) d J k 2 + x~ (4n)^^ 2' 

and we define d = 4 — 2e. (We work in Euclidean space throughout.) The evaluation of 
I(x, y, z) is non-trivial, and will be the subject of Section 4. 

We turn first to the subject of renormalisation. We choose to work throughout with 
renormalised parameters gz, g, g', m 2 , h, and A. ((73, g, g' are the three gauge couplings, 
with the usual conventions). Then rather than compute separately a set of one-loop dia- 
grams with counter-term insertions, we subtract (minimally) the sub- divergences diagram 
by diagram . is then obtained by simply taking the finite parts of the resulting ex- 

pressions, discarding the 1/e 2 , 1/e poles which are automatically cancelled by the usual 
renormalisation constants (which we need not calculate). For graphs not involving vector 
bosons, this procedure amounts simply to the replacement of I(x, y, z) and J(x, y) by their 
subtracted values: 

I(x, y, z) -> I(x, y, z) = I(x, y,z)-— (J(x) + J{y) + J{z)) (3.4) 

J(x,y) -> J(x,y) = J(x,y) + —(xJ(y) +yJ(x)). (3.5) 

Ke 

There is one complication, however. When vector bosons are present, the algebra involved 
in reducing the graph to dependence on / and J may produce explicit factors of d. This 
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gives rise to an apparent ambiguity: what is the subtracted form of dl? The answer is that 
the result depends on which subgraph produced the factor of d. Thus when vector fields 
are present we must explicitly evaluate the contribution of the subtractions to each graph. 
This is still, however, much easier than explicitly considering counter-term insertions. 
Because we work in the 'tHooft-Landau gauge, the gauge parameter is unrenormalised 
and so creates no difficulties. 

We conclude this section with the results for J which will be relevant subsequently. 
It is straightforward to show from (3.3) that: 

k 2 J(x, y) = — + xy (l — lna; — lay + lmclny) , (3.6) 
K J(x) = ~+x(Ex-l), (3.7) 
K 2 eJ(x, y) = — + xy (2 — Inx — Iny) . (3.8) 

It is the finite part of each expression that is substituted in our expressions for in 
Section 5. Terms of the form of Eq. (3.7) and Eq. (3.8) appear in connection with the d 
dependence discussed above. 



4. Evaluation of I(x,y,z) 

In this section we indicate how the differential equations method in conjunction with 
the method of characteristics leads to a simple form for the integral I(x, y, z). Of course 
there has been much effort expended on higher loop Feynman integrals, and many powerful 
techniques have been developed (for some examples see ref. [14]); nevertheless we feel that 
our method here has some interesting features and is worth presenting in detail. 

We start from the identity: 

/d 
d d kd d 
dk^ (k 2 + x)(q 2 + y)((q + k) 2 + z) 

From Eq. (4.1) it follows that 

2x— + (z + x-y)— = (d-3)I + K x {x, y, z) (4.2) 

where 

K 1 (x,y,z) = -^-(J(x)-J(y)). 

Now from Eq. (4.2) and similar equations produced by (x,y,z) permutations, one can 
eliminate dl/dy and dl/dz and then solve the resulting equation in a similar way to that 
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adopted in ref. [7]. It turns out, however, that a more elegant solution follows if we start 
with the following equation: 



where 



. .dl . ,9/ , . dl _ , 



K(x,y,z) = Kxix.y.z) +K 1 (y,z,x) + K 1 (z,x,y) 



= -((z-x)(zx) e + (x-y)(xy) e + (y - z){yz) e )V 



(4.3) 



(4.4) 



and 



d. 



r' = (^) 2 T(2-^)r(i-^)(47r)- 

The method of characteristics involves solving a system of ordinary differential equa- 
tions, as follows: 

dy 



dt 



dx 



dz dl 



y — z z — x x — y K 



(4.5) 



subject to initial conditions which we shall choose to be x = X, y = Y and z = (at 
t = 0), and we will suppose without loss of generality that X > y.We then have: 



I(x, y, z) = I(X, Y, 0) + f dt' K(x(t'),y(t'),z(t')). 

Jo 

Using Eq. (4.4) and Eq. (4.5), we can rewrite Eq. (4.6) as 



(4.6) 



I(x,y,z) = I(X,Y,0)-T' 



dx (yz) 6 + / dy (zx) 6 + / dz (xy)~ 
x Jy Jo 



(4.7) 



Now from Eq. (4.5) it follows that for all t 



x 2 + y 2 + z 2 = d 2 = X 2 + Y 2 
x +y + z = c= X +Y 



(4.8) 
(4.9) 



where c and d are constants. Therefore 



xy = z —cz+-(c —d) 
with similar equations for yz and zx. Hence Eq. (4.7) becomes 



I(x,y,z) = I(X,Y,0)-V 



2 f2 

-y 



ds(s 2 -a 2 )- £ 



(4.10) 



(4.11) 
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where 



a 



y - j = \{X - Y) = l -{x 2 +y 2 + z 2 - 2xy - 2yz - 2zx)^ . 



(4.12) 



Now I(X, Y, 0) is also tricky to evaluate by elementary methods; but by employing the 
method of characteristics once again, one can show that 



i(x : y, o) = i(x - y, o, o) + r' / ds \s 2 --(x- y) 2 ] \ 

J\{X-Y) ' 4 



(4.13) 



Substituting Eq. (4.13) in Eq. (4.11) and using Eq. (4.9) and Eq. (4.12) we obtain 



/(,•. u, z) = I (2a. 0, 0) + I'' F(°- -y) + F(| - z) - F(x - °-) 



(4.14) 



where 



F(w) = / ds(s 2 -a 2 )- e . 

J a 



(4.15) 



Since I(2a, 0, 0) can be evaluated by elementary methods, we have reduced the problem to 
a single integral, F(w). In spite of appearances, Eq. (4.14) is symmetric with respect to 
(x, y, z) permutations, since it is easy to show that, for example, 

F(l-y)-F(x- C -) = F( C --x)-F(y- C -). 

However Eq. (4.14) is only valid in the region a 2 > 0. In the region a 2 < 0, it is possible 
to derive the following form of the solution: 



I(x, y, z) = -I(2b, 0, 0) sin ^ + V [g(°- - x) + G(°- - y) + G(°- - z) 



(4.16) 



where 



G(w)= ds(s 2 + b 2 )- £ and b 2 = -a< 
Jo 



(4.17) 



It is a nice exercise to show that for z = x, Eq. (4.16) can be rewritten as Eq. (11a) of 
ref. [7]. Note that for a 2 = 0, which in x,y,z space is a cone with its apex at the origin, 
the integral is trivial. 
Writing 

(s 2 + b 2 )~ e = 1 - e ln(s 2 + b 2 ) + \e 2 \n 2 (s 2 + b 2 ) + ... (4.18) 



and using 



(4^,0,0) ■ £gzigyazg ( - r 



(4.19) 
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we obtain from Eq. (4.16) that 



k 2 I{x, y, z) = - ^ - ^ (y - Li) 



- - 1/2 - 6L1 + (y + z - a^lnylnz 
+ (z + x — y)\nz\nx + (y + x — z)lny\nx 
+ ^(a:,y,z) + C (7 + C(2)) 



where 



- — m - — m 



L m = xln x + y\n y + zln z 



and 



y, z) = 86 [L(^) + L(^) + L(9 Z ) - | In 2] . 
L(£) is Lobachevskiy's function[15], defined as 



L(t) = — dx In cos x. 
Jo 



(4.20) 

(4.21) 
(4.22) 

(4.23) 



The angles 9 X , 6 y , 9 Z are given by 



9 X = tan 1 ( 2 X ) etc. 
b 



(4.24) 



Eq (4.22) is valid only in the region a 2 < (ie. inside the cone). For a 2 > 0, we 
obtain from Eq. (4.14) a result identical to Eq. (4.20) except that now 



y, z) = 8a \-M{-cj> x ) + M{(j> y ) + M(0,)] 



(4.25) 



where 



M(t) = - I dx lnsinh 
Jo 



x 



and 



are given by 



coth 1 (1 -) etc 



a 



£(x, y, z) is /U-independent and therefore plays no role in the RG analysis of section 6. 

It is interesting to compare the results Eq. (4.20) and Eq. (4.22) with those of ref. [7] 
for the special case z = x. In that case we have 

sin# x = sin# z = ]-{-) 5 



sin B v = 1 

y 2x 



and using the identity 

2L^sin -1 (^y)^ +L^sin- 1 (^l - = |ln2 + 2^ ' * ^ ln(2 sinx) dx (4.26) 

it is easy to show that Eq. (4.20) reduces to Eq. (12) of ref. [7]. A similar process works 
in the a 2 > case. 

We will make extensive use of the form / takes when its subdivergences are subtracted, 
I. From Eq. (3.4) we find that 

/ = / + ^2( c +( c - L e +{\l*-Li + c+ ^<(2))e 2 + ...). (4.27) 
It is the finite part of this expression which we use subsequently. 



5. Two-loop results. 

In this section we present our results for the various contributions to in (what we 
hope is) a clear and systematic manner. It is natural to divide the calculation into parts 
according to the nature of the contributing fields, thus 

V<® =V S + V SF + V sv + V FV + V v (5.1) 

where S, F, V denote scalar, fermion and vector fields respectively. 

The Feynman rules of the standard model are well known; for the calculation of the 
effective potential we must only remember that we are not calculating at the tree minimum 
of the potential. Apart from giving the G, bosons a non-zero mass, this makes little 
difference. The results are as follows: 

V S =Z ^f- [HH, H, H) + I(H, G, G)] 

+ ^[J(H,H)+2J(H,G) + 5J(G,G)} (5.2) 

o 

(in fact this is the special case = 4 of the calculation reported in ref. [7]). 

V SF = 3/i 2 {(2T - ^H)I(T, T, H) - \gI(T, T, G) + (T - G)I(T, G, 0) 

+ J(T, T) - J(T, H) - 2 J(T, G)}, (5.3) 
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V sv = TT^T-nMH, G, Z) + g2{1 2 ™ 9)2 A(G, G, Z) + \e 2 A{G, G, 0) 
8 cos 2 8 cos 2 2 

+ -g 2 (A(H, G, W) + A(G, G, W)) - g 2 sin 4 6ZB(Z, W, G) 

- e 2 WB(W, 0, G) - \g 2 WB(W, W, H) - ^-^ZB{Z, Z, H) 

g 2 (l-2sm 2 9) 2 c f_ 

4 cos 2 (9 8cos 2 6> v v " 

+ l g *(C(W,H) + 3C(W,G)), (5.4) 
a 2 

V v = - — {2 sin 2 0A(W, W, 0) + 2 cos 2 6A(W, W, Z) - 2£(W, 

- 4 cos 2 0E(W, Z) + 2A(0, 0, + cos 2 0A(O, 0, Z)} (5.5) 
(The last two terms represent the ghost contribution). 

V FV = -3j2 M + a )) D ^ Z) + (v 2 - a})E(F, F, Z)] 
f 

- ^g 2 (n G - 1)15(0, 0, W) - \g 2 n G D{^ 0, W) - \g 2 D(T, 0, W) 

- {Ag 2 + ^e 2 ) [D(T, T, 0) + E(T, T, 0)] . (5.6) 

In Eq. (5.6) the sum over / is over all quarks and leptons, and Vf and a/ denote the vector 
and axial couplings to the Z boson. Thus, for example, 

ff(l-fsin 2 fl) g 

v t = and a t = - 5.7 

4 cos t) 4 cos v 

where 6 is the usual weak mixing angle. F = except for the top quark, when F = T. 
no is the number of generations; of course since we neglect all Yukawa couplings except h, 
our calculations are really only applicable for n G = 3. The various functions introduced 
in Eqs. (5.4)-(5.6) are defined as follows: 

A(x,y,z) = -{-4a 2 I(x,y,z) + (x - y) 2 I(x,y,0) + (y - x - z)J(x,z) 

z 

+ (x-y- z)J(y, z) + zJ(x, y) + 2z(x + y- •/(*)}, (5.8) 

B(x, y, z) = -— {(lOxy + z 2 + x 2 + y 2 - 2xz - 2yz)I(x, y, z) 
4xy 

+ (2zx - x 2 - z 2 )I(x, z, 0) + (2zy - y 2 - z 2 )I(y, z, 0) 

+ z 2 I(z, 0,0) + (x + y- z)J(x, y) 

-yj(x,z) -xJ(y,z) + 6xyK~ 1 (J(x) + J(y)) 

-8xyel(x,y,z)}, (5.9) 
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C(x, y) = 3J(x, y) - 2eJ(x, y) - 2yK~ 1 J{x), (5.10) 

D(x, y, z) = -{-(x 2 + y 2 - 2z 2 +xz + yz- 2xy)I(x, y, z) + (x - y) 2 I(x, y, 0) 

z 

— 2zJ(x, y) + (2z + y — x)J(x, z) 

2 

+ (2z + x- y)J{y, z) + -z(2z - 3y - 3x)k~ 1 J(z) 

+ 2ez[(x + y- z)I(x,y,z) + J(x,y) - J(y,z) - J(x,z)}, (5.11) 

E(x, y, z) = ^(6I(x, y, z) + 4k" 1 J(z) - 4el(x, y, z)), (5.12) 
A(x, y, z) = A(x, y, z) + A(y, z, x) + A(z, x, y) (5.13) 



where 



A(x, y, z) = {(x 4 — 8x 2 yz — 2x 2 y 2 + 32a 2 xy)I(x, y, z) 

-((x 2 - y 2 ) 2 + 8(x - y) 2 xy)I(x, y, 0) + x 4 I(x, 0, 0) 

+ zJ(x, y)(9x 2 + 9y 2 — 9xz — 9yz + 13xy — z 2 ) — (— — l)xyzJ(x 7 y) 

+ 8(x — y) 2 xyel(x 1 y : 0) + 8z(xz + yz — x 2 — y 2 — xy)eJ(x, y) 
25 

— Axyz{ — x + 6y + Qz)k~ 1 J(x) — 32a 2 xyel(x, y, z)}, (5-14) 
o 



and 

E(x,y) = ^J(x,y) + " ™)j(x,y) - \n~\xJiy) + yJ(x)). (5.15) 

The a 2 that appears in Eq. (5.8), Eq. (5.14) is defined in Eq. (4.12). The apparent singu- 
larities in, for example, A(x, y, z) for z = are easily seen to be spurious. 

It is straightforward (but tedious) to substitute Eqs. (5.8)-(5.15) in Eqs. (5.2)-(5.6). 
The resulting expression is not particularly illuminating, however, and so we do not present 
it. Clearly evaluation by computer will be necessary when we come to applications, in any 
event. 



6. The renormalisation group. 

V{4>) obeys the following RG equation: 
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where Aj = {m 2 , A, g, g', g$, h}. We begin by verifying Eq. (6.1) at leading order. We have 



p(i)y(0) = ^lyti) ( 6 .2) 



where 

aA, ' 

V 

Using Eq. (2.4), Eq. (6.2) becomes: 



kV W y(0) = hff 2 + 3G 2 - 12T 2 + 6W 2 + 3Z 2 ) . (6.3) 

Note that there is a (^-independent term on the RHS of Eq. (6.3). We therefore need to 
add a suitable term to of the form f(Xi)m 4 . Equivalently, of course, one can simply 
redefine by subtracting it at = 0[16]. This just shifts the potential by a constant, 
and has no consequences for the considerations of this section. It will, however, contribute 
non-trivially to the RG-improved form of the potential [17]. 
Comparing coefficients of (j) 4 and </> 2 , we obtain 

Av(/^ 1} - 4A 7 (1) ) = ^(16A 2 - 144/i 4 + 9g' 4 + 18g 2 g' 2 + 27g 4 ) (6.4) 

and 

k(/3 {1) 2 - 2m 2 7 (1) ) = 2m 2 A. (6.5) 



Substituting 



^ = \(12h 2 -9g 2 -3g' 2 ) (6.6) 



we obtain the well known results 

K pW = 4A 2 + 12Xh 2 - 36/i 4 - 9Xg 2 - 3Xg' 2 

+ -/ 4 + \g 2 g' 2 + ^g 4 (6.7) 

and 

K^l=m 2 (2\ + 6h 2 - 9 -g 2 - 3 -g' 2 ). (6.8) 
At the two-loop level we have 

p(2)y(0) = _ ( 6 . 9 ) 
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The evaluation of (fj,d/dfi)V( 2 ) using the results of section 5 is a straightforward application 
of the following results: 





= 2(Li - 


•2c), 


(6.10a) 


k 2 ^-^-J(x, y) 
d/i 


= 2xy(2 - 


- hue — lmy), 


(6.106) 


d 


= -2c, 




(6.10c) 


d 


= 4xy, 




(6.10d) 


K^ — Jyx) 


= -2x. 




(6.10e) 



Using Eq. (6.10) and Eqs. (5.2)-(5.6) we find that 

K 2 /i|V 2) = <r 4 {^A 3 + 2ft 2 A 2 - 7ft 4 A - 12ft 6 + 5ft"</ 2 + 32 9 2 ft 4 + 9 ' 4 ft 2 

" 3 ^ + |V + ^ V + 1a 9 2 9 « - fA 2 , 2 - iA 2 ,' 2 - f J 

+ ^ + I 9 v 2 + 5 9 y 4 + 5 9 <») + § 9 y 2 + if 9 y 4 + g n 

+ m 2 2 {4A 2 + 12A/* 2 - 18/i 4 - ^A(3# 2 + g' 2 ) 

~ Y6 g4 + Y g2g ' 2 ~ T6 9 ' 4 ^ + lo S arithmic terms. (6.11) 

The logarithmic terms in Eq. (6.11) are not given explicitly because they simply 
cancel an identical set of terms from V^V^K This provides an excellent check on the 
calculation, since there are terms of the form In if, InG, InT, InW and InZ which must 
cancel separately. 

The calculation of p( 1 )y( 1 ) is straightforward; the necessary one-loop RG functions 
are given by Eq. (6.6)-(6.8) and: 

«/V = \^i\ n G + = 9*{\n G - ^). (6.12) 

We hence obtain from Eq. (6.9) that 



k 2 ((3 {2) - 4A 7 ^) = --A 3 - 24A 2 /i 2 + 6A 2 (3# 2 + g' 2 ) + 24A/* 4 + ^Ao 4 



+ ^Xg 2 g' 2 + |a ' 4 + ISO/* 6 - 192/^f - 16h*g» - y/>V 

+ 63 W - f*V 4 + 3 «x " 8710)96 " ( 1 + f "^V 2 

,239 40 , o /4 ,59 40 , , 6 , 
" (-24 + y^)^V 4 - + Y n G )g' G } (6.13) 
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and 



k 2 (/3^ -2mV 2) ) =2m 2 {-A 2 -6Ah 2 + 2A(3 9 2 + 9 ' 2 ) 



+ l^ 4 + !^' 2 + f|^}. (6.14) 

We now require r )^ 2 \ For a general gauge theory (and in a general covariant gauge) 
this is given in ref. [18]; we have merely to specialise to the Landau gauge and the SM, 
where care must be taken to keep track of the factors. The result is: 

K 2 (2) * A 2 _ 27^4 20 2 h 2 45 2 2 85 ^ 

'6 4 3 8 24 

+ (f"°-f)« 4 + ^ + (l,o + |)^ (6.15) 

Substituting Eq. (6.15) in Eq. (6.13) and (6.14) we obtain the final results: 

2fl( 2 ) - _ 0/1X27,2 , fi \2/o 2 , „/2n _ o U 4 

_i_ sn\^2i,2 



= - - 24A^ + 6A W + ^) - 3A/i 4 + 80A^/i 

+ | Xg 2 h 2 + ^-Xg' 2 h 2 + (lOn G - ^)A<? 4 + ^V</ 2 



+ (yn G + ^)A</ 4 + ISO/. 6 - 192h 4 g 2 - 16/. V 2 - y*V 

+ 63/W - y/^V 4 + 3((^ - 8n G )g« - (g + j^Y 2 
,239 40 x , ,4 /59 40 wfi \ , 

- ^ + y v " ( 21 + t no)g ) • (fU6) 



and 



k 2 /^ = 2m 2 {-^A 2 - 6A/. 2 + 2A(3^ + g' 2 ) - ^h 4 + 20g 2 h 2 + fg 2 h 2 

85 /2j2 /5 385\ 4 15 2 , 2 ^25 157\ , 
+ + ( 5 « G - ^ j / + ^V 2 + + — j 9 «}. (6.17) 

We can compare our result for /^ 2) , Eq. (6.16), directly with Eq. (B.8) of ref. [6]. Agreement 
is complete apart from 

(i) the sign of the Xg 2 g' 2 , Xg' 4 terms. 

(ii) the magnitude of the Xg' 4 term. 

From the result of a general gauge theory (Eq. (4.3) of ref. [6], or see also ref. [19]) it 
appears that these discrepancies arise from errors in the reduction to the SM case in ref. [6] 
rather than in the general result. f Unfortunately the incorrect form of Eq. (B.8) has been 
applied by other authors to running coupling analyses, e.g. ref. [20]. The numerical effect 

(2) 

of the error on these calculations is probably small, however. The result for j3 2 can be 



| See note added. 
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verified from the general formulae of ref. [19] and will be needed when we consider the RG 
improved form of V. 

At this point it is appropriate to comment on the work of Alhendi[5]. He essentially 
reverses the procedure adopted in this section, in order to deduce the //-dependent terms 
in given the MS /9-functions. Unfortunately his expressions for 7^) and (3^ 2 are 
incorrect; he assumes for example (in our notation) that the coefficient of the g 4 term in 
7^ 2 ) is just \ of the coefficient of the Xg 4 term in (3^\ This amounts to the assertion 

that there are no 1PI contributions of this kind to (3^\ which is not correct (even in the 
Landau gauge). Aside from this comparison would still be difficult since he does not use 
the MS form of V^. These problems aside, it is clear that by this method one can derive 
the form of the potential in the region (p 2 ^> |m 2 | only. 



7. Conclusions and outlook 

We have presented a calculation of (4>) in the standard model, using dimensional 
regularisation and minimal subtraction. Applying the renormalisation group to the result 

(2) (2) 

led to numerous checks and also expressions for j3\ ' and .We used differential equations 
and the method of characteristics to find the relevant Feynman integral; an approach which 
seems to us preferable to more traditional techniques. 

Studies of the stability of the electroweak vacuum reviewed in ref. [2] suggest a limit 
of around m t < 95GeV + 0.6m# for the top and Higgs masses. Our calculation will enable 
us to probe further the sensitivity of this result to radiative corrections. We will report 
these calculations elsewhere; the following argument, however, suggests that it will not be 
affected dramatically: Instead of solving the RG equations to produce an "improved" V, 
suppose we take the unimproved V and simply set fx = 4>. All couplings (and m 2 ) then 
become functions of 0; but this choice of \i means that when 4> 2 ^> |m 2 |, the radiative 
corrections to V, although non-zero, have no large logarithms. As long as the running 
couplings remain perturbative, we may therefore expect the radiative corrections to cause 
only a small change in the above limit. Perturbative confidence in the limit means that 
should it turn out that m t and m# fail to satisfy it, this would provide compelling evidence 
for physics beyond the standard model. 

Note Added 

Our thanks to Mike Vaughn for correspondence, and for confirming that the authors 

(2) 

of ref. [6] now agree with our result for (3\ \ Eq. (6.16). 

The published version of this paper (and unfortunately also the published erratum) 
contain a number of typographical errors and omissions, which have been corrected here. 
These errors affected the result for (Eq. (5.1)) and the result for (3^1 (Eq. (6.17)), 
but not the result for /^ 2) (Eq. (6.16)). 
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We are grateful to Peter Arnold, Steve Martin and Mingxing Luo who drew most of 
these errors to our attention. See ref. [21]-[23]; ref. [22] includes a generalisation of the two 
loop effective potential calculation to the case of an arbitrary gauge theory. 
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